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Abstract 


This  report  deals  with  the  propagation  of  a  pulse 
in  a  semi-infinite  inhomogeneous  medium  for  which  the  law 
of  variation  of  the  refractive  index  is  such  that  the  wave 
equation  can  be  solved  by  the  method  of  separation  of  vari- 
ables. The  Laplace  transform  of  the  solution  is  first  ex- 
pressed as  a  series  of  eigenf unctions,  whose  asymptotic 
behavior  at  distant  points  of  the  Laplace-plane  is  inves- 
tigated by  means  of  improved  variants  of  the  W.K.B. 
method.  An  approximation  is  found  whose  inverse  Laplace 
transform  can  be  calculated}  it  applies  when  a  shadow,  in 
the  usual  sense  of  geometrical  optics,  is  formed,  and  it 
describes  the  variation  of  the  solution  immediately  behind 
the  diffracted  front  which  propagates  into  the  shadow. 

The  work  is  an  extension  of  the  author's  recent  report  on 

\3\ 
diffraction  of  pulses  by  a  circular  cylinder  L  J . 
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Introduction 

In  a  recent  report'-'3-'   the  author  discussed  the  diffraction  of  two- 
dimensional  pulses  by  a  circular  cylinder.     The  object  of  the  present  report 
is  to  extend  this  discussion  to  a  more  general  class  of  propagation  and  scat- 
tering problems. 

Let  0  be  a  function  of  two  space  coordinates  x  and  y  and  of  the 
time  t,  which  satisfies  the  two-dimensional  wave  equation 


d20 
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c     at 

where  c,  the  velocity  of  wave  propagation,   may  be  a  function  of  x  and  y.     This 
is  the  equation  satisfied  by  the  pressure   (or  its  time  integral)  of  an  acoustic 
disturbance  in  an  atmosphere  whose  temperature  is  not  necessarily  uniform,    or, 
if  c  is  constant,   by  the  transverse  components  of  TE  -  or  TM-waves.     One  may 
consider  the  propagation  of  waves  governed  by  this  equation  in  the  space  exterior 
to  a  cylinder  with  generators  parallel  to  the  z-axis,   on  whose     surface  either 
0  or  its  normal  derivative     is  to  vanish.     But  it  is  always  possible  to  effect 
a  conformal  transformation  of  x  and  y  which  maps  this   'physical  space'   into  the 
upper  half-plane,   y  >  0,   while  preserving  both  the  form  of  the  wave  equation 
and  of  the  boundary  condition.     It  is  therefore  sufficient  to  consider  the  pro- 
pagation problea  in  an  inhomogeneous,   semi- infinite  medium  with  either  of  the 
boundary  conditions  0  -  0  or  80/dy  -  0  at  y  -  0.     The  basic  problem  is  the 
determination  of  the  Green's  function,   G,   which  can  be  defined  as  the  dis- 
turbance generated  by  an  instantaneous  point  source.     Using  Dirac's  5-f unction, 
one  can  say  that  the  Green's  function  satisfies  the  equation 
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6(x  -  xo)6(y  -  yQ)6(t),         (yQ  >  0) 


and  also  the  initial  conditions 


G  -  0,  H  -  0,   when  t  -  0,  (x  -  xq)2  +  (y  -  yQ)2  /  0  , 

with  the  appropriate  boundary  condition  at  y  ■  0.  G  is  the  disturbance  caused 
by  an  instantaneous  source  at  x  =  x  ,  y  =  y  at  the  time  t  «*  Oj  the  disturbance 
due  to  such  a  source  acting  at  time  t  is  of  course  simply  G(x,y,t  -  t  ,x  ,y  ) . 

If  the  medium  were  of  infinite  extent,  the  Green's  function  would 
reduce  to  Hadamard's  elementary  solution  divided  by  2n.  Let  p  =  p(x,yjx  ,y  ) 
be  the  time,  as  determined  by  Jermat's  principle,  when  a  wave  starting  at 
t  =  0  from  (x  ,y  )  arrives  at  the  point  (x,y) .  Then  the  elementary  solution, 
G  ,  vanishes  for  t  <  p,  and  is  given  at  least  for  small  positive  t  -  p  by  a 
series  of  the  form 

G0  ■  22  Vx>y)  (t  -  p  >        • 

°    n=0  n 

The  coefficient  E  can  be  determined  recursively,  essentially  by  integration 
of  the  wave  equation  along  the  ray  from  (x  ,y  )  to  (x,y)  . 

If  there  is  a  boundary  at  y  =  0  then  a  ray  from  (x  ,y  )  meets  y  =  0 
at,  say,  (x^  ,0)  ,  and  one  can  construct  a  ray  reflected  from  (x..,0)»  If  (x,y) 
is  a  point  on  this  reflected  ray,  then  one  can  determine  the  time  required 
for  a  wave  to  travel  from  (x  ,y  )  to  (x^ ,0)  along  the  direct  ray,  and  from 
(x,  ,0)  to  (x,y)  along  the  reflected  ray.  This  total  time  is  a  function  of 
position,  and  will  be  denoted  by  p  .  Now  it  is  possible  to  construct  an- 
other  solution,  G  ,  of  the  wave  equation.  This  solution  vanishes  for  t  <p  , 
and  for  t  >  p  it  is  of  the  form 
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*    J2°   *        2   *2  ""2 

G  ■  2H  En  (x»y)   (*  -  P  )    5 

n=0 

further,  G  is  such  that  on  y  =  0  either  G  =  -  G  or  dG  /dy  =  -  dG  /dy.  The 

•s- 
E  are  obtained  recursively  by  integration  of  the  wave  equation  along  a  re- 
flected ray.  Clearly, 

G  =  G  +  G* 
o 

satisfies  the  conditions  of  the  problem.  The  leading  coefficients,  E  and 
E  ,  can  be  obtained  easily;  the  approximation  which  results  when  only  the 
leading  terms  are  retained  corresponds  to  the  results  obtained  by  the 
methods  of  geometrical  optics. 

However,  it  may  happen  that  there  are  points  of  the  xy-plane  which 
cannot  be  joined  to  (x  ,y  )  by  means  of  a  direct  and  a  reflected  ray.  All 
the  domains  whose  points  have  this  property  are  called  the  shadow,  and  in 
this  case  we  shall  say  there  is  diffraction.  In  the  shadow  neither  p  nor 
the  E  can  be  defined,  and  there  arises  the  problem  of  finding  an  alternative 
representation  of  G. 

The  present  report  gives  a  partial  solution  of  this  problem.  It 
is  confined  to  the  case  where  the  refractive  index  is  of  the  form 


\   ■  *(*)  +  g(y)  , 


c 

where  the  functions  f (x) ,  g(y)  are  subject  to  certain  restrictions,  the  most 
important  of  which  is  that  g(y)  must  be  increasing  and  tend  to  infinity  with 
y.  An  example  discussed  in  Appendix  II  suggests  furthermore  that  the  solution 

obtained  may  not  be  valid  in  the  whole  shadow. 

f3l 
The  method  used  is  similar  to  that  used  in  the  previous  report L  J . 

First,  the  Laplace  transform  0f  of  <j(,   namely 


-u  - 


0    =   /    e"st  0  dt  , 


is  expanded  as  a  series  of  eigenf unctions  which  represent  modes  propagated 
parallel  to  the  boundary  (Section  3) •  Next,  the  asymptotic  behavior  of  the 
eigenvalues  as  Re(s)  -^  oo  is  investigated.  This  is  done  by  a  variant  of 
Langer's  method  (Sections  h   -  6) ,  which  was  used  in  the  cylinder  diffraction 
problem.   [An  alternative  method  has  recently  been  found  and  is  described 
briefly  in  Appendix  I.  This  method  is  simpler  and  -  at  least  from  a  formal 
point  of  view  -  more  effective.]  It  is  then  possible  to  derive  approxima- 
tions valid  for  Re(s)  large  for  the  terms  of  the  series  for  0*     The  form  of 
these  approximations  suggests  that  they  are  valid  in  the  shadow  only,  and 
shows  that  only  the  leading  term  of  the  first  mode  may  be  retained  (Sections 
7  and  8).  In  Section  9  the  results  obtained  are  related  to  the  geometrical 
optics  of  the  problem.  The  final  step,  namely  the  calculation  of  the  in- 
verse Laplace  transform  of  the  approximation  obtained  for  0,  is  effected  by 
appealing  to  results  already  obtained  in  the  case  of  diffraction  by  a  cir- 
cular cylinder.  Since  that  case  is  typical  of  the  case  of  diffraction 
problems  we  are  considering>  the  present  report  repeats  the  methods  and 
calculations  used  in  the  earlier  work;  however,  the  results  obtained  are 
more  general  and  can  be  used,  for  instance,  to  calculate  the  diffracted 
pulse  in  the  shadow  of  an  elliptic  or  a  parabolic  cylinder  (see  Section  2). 
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1.  -formulation  of  the  problem 

The  problem  can  be  formulated  analytically  as  follows.  It  is  re- 
quired to  find  a  function  0  of  two  variables  x,y  (the  space  coordinates)  and 
the  time  t  which  satisfies,  in  the  domain  y  >  0,  the  wave  equation 

2  2  2 

(1.1)  *-$     +     *!$=     [f(x)   +  g(y)]  §-f    -     6(x  -  XQ)6iy  -  yQ)  6(t),      y>0 
dx  dy  dt 

and  the  initial  conditions 

(1.2)  0     =  |?  =  0,    t  =  0,  y  >  0,  (x  -  x  )2+  (y  -  y  )2  /  0  . 

In  addition,   either  of  the  boundary  conditions 

(1.3)  M       -     o 

LJy=0 

or 

is  to  hold  at  y  =  0.  The  symbol  6  in  (1.1)  stands  for  the  Dirac  6-function. 
The  functions  f(x)  and  g(y)  are  assumed  to  satisfy  such  conditions  of  reg- 
ularity (or  analyticity)  as  the  argument  may  require,  as  well  as  the  follow- 
ing conditions: 

(1.5)  g(0)  >  0  ; 

(1.6)  f(x)  >  0  ; 

(1.7)  g  (y)  >  0  for  y  >  0  j  g(y)  -}  oo  as  y  -*■  co  . 

2c  Illustrative  examples 

As  already  mentioned  in  the  Introduction,  problems  of  the  type 
just  described  arise  in  the  discussion  of  the  diffraction  of  two-dimensional 
pulses.  To  illustrate  this,  let  us  consider  the  following  situation  for 
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definiteness.  Suppose  that  0  is  the  velocity  potential  of  an  acoustic  dis- 
turbance. Consider  first  a  homogeneous  medium  in  which  the  space  coordinates 

i   i 
are  (x  ,y  )  ;  then  we  have 

(2.i)  i^  ♦   £|    .  jj  a?g   , 

3x       3y       c   dt 

where  c  is  a  constant.  We  may  then  consider  the  field  due  to  a  pulse  re- 
presented by  Hadamard's  elementary  solution 

1    H(ct  -  R)       2      i   i  2     «    i  2 
(2.2)      j*0=  2J   2  2       Vg  ,   R  -  (x-x)  ♦  (y-y') 

(c  t  -  R)  ' 

(H  denotes  Heaviside's  unit  function).  This  field  may  be  derived  by  including 
a  6-f unction  'source'  term  in  (2.1)  analogous  to  that  in  (l.l),  and  adding 
initial  conditions  analogous  to  (1.2).  Let  the  reflector  (scattering  object) 
be  bounded  by  a  surface  S  at  which  either  0  ■  0  or  d$/dn   =  0,  where  d/dn  de- 
notes differentiation  along  the  normal.  It  may  happen  that  a  conformal  map- 
ping of  the  exterior  of  S  onto  the  upper  half  of  the  xy-plane  transforms 
(2.1)  into  an  equation  of  the  type  (1.1);  then  the  scattering  problem  is 
equivalent  to  that  formulated  in  Section  1. 

Example  I.  Circular  cylinder.  The  conformal  mapping 

'   .  '     -i(x  +  iy) 
x  +  iy  =  e       J 

•2   '2 

transforms  the  exterior  of  the  unit  circle  x  +  y  =  1  into  y  >  0,  and  (2.1) 

into  (1.1)  with  f (x)  =  0,  g(y)  =  c~  e  y.  In  this  case,  0  must  be  a  periodic 
function  of  x  with  period  2n.  But  it  can  be  expressed  as  the  sum  of  terms, 
each  of  which  satisfies  the  conditions  of  the  problem  of  Section  1. 
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Example  II*  Parabolic  cylinder*  The  conforraal  mapping 
X  +  iy  =  ^  (x  +  ly) 

'2      • 

transforms  the  exterior  of  the  parabola  y   >  U(x  +  l) ,  into  the  region  y  >  s/2. 

2  2         2  2 
At  the  same  time,  (2.1)  becomes  (l.l)  with  f(x)  =  x  /c  ,  g(y)  =  y  /c  .  The 

reduction  of  the  scattering  problem  to  that  of  Section  1  -  by  a  translation 

of  the  origin  along  the  y-axis  and  a  change  of  scale  -  is  trivial. 

Example  III.  Diffraction  by  an  elliptic  cylinder  .  Let  the  mapping  be 

t    i 

x  +  iy  =  sin(x  +  iy)   • 

—2  —2 

Then  f (x)  »  c   cos  2x  and  g(y)  =  c   cosh  2y,  and  y  >  a  corresponds  to 

i2  i2 

2~  +    — 2—  -   l  ♦ 

cosh  a      sinh  a 

This  again  gives  essentially  the  problem  of  Section  1,  with  the  same  additional 
difficulty  that  occurs  in  the  case  of  Example  I  (circular  cylinder);  this  can 
be  circumvented  in  the  same  manner. 
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3.     The  formal  solution 

To  obtain  a  formal  solution  of  the  problem  stated  in  Section  1,  we 
begin  by  introducing  the  Laplace  transform  0  of  0,  i.e., 

(3.1)  ?(x,y,s)      =     /     e"st     0(x,y,t)dt  . 
This  satisfies  the  equation  (see    [3]   p.  6,  footnote) 

(3.2)  2-f     ♦     2l|      =       s2[f(x)    ♦  g(y)l   ?     -     6(x-x)6(y-y  ) 

dx  dy  J  00 

and  either  of  the  boundary  conditions  0  =  0  or  30/dy  =  0  at  y  =  0.  In  addition, 
as  a  consequence  of  the  general  properties  of  Laplace  transforms,  0  -^  0  as 
Re(s)  -^oo.  fin  practice  this  appears  to  be  equivalent  to  the  'radiation  con- 
dition  0->Oasy-^oq/.  Since  this  last  condition  may  be  taken  to  be  a 
consequence  of  the  assumption  that  0  =  0   for  t  <  0,  we  shall  refer  to  it  as 
the  'pulse  condition'. 

We  can  now  obtain  a  formal  solution  of  (3.2)  as  an  eigenf unction 
expansion  by  using  the  method  of  normal  coordinates  of  the  theory  of  vibra- 
tions. Consider  the  solutions  of  the  ordinary  differential  equation 

2 

(3.3)  ^  +  (X  -  s2  g(y))  X  -  0 

dy 

where  X  is  a  real  or  complex  parameter.  For  sufficiently  large  y,  the  coef- 
ficient of  %   is  negative,  and  there  are  two  independent  solutions  cf  (3«3) 
which,  for  large  s,  are  asymptotic  to 

[s2g(y)  -  x]      exp  J  +  /  (s2g(y)  -  X)l/2  dy 


If  we  require  %   — ^  0  as  Re(s)  —^co,  then  we  must  choose  the  minus  sign,  and 
the  solution  thus  determined  has  the  property  that  it  tends  to  zero  as  7  «^co. 
LetX(y»s>^)  denote  such  a  solution  of  (3«3)  .  Since  g(y)  ->oo  as  y  -}  00,  this 
solution  is  determined  except  for  a  factor  independent  of  y«  Each  of  the 
boundary  conditions  X  ■  0  at  y  =  0  or  X  =  0  at  y  =  0  is  satisfied  only 
for  a  sequence  of  values  of  X,  say  X.  ,  X~  ...,  the  eigenvalues  of  the  problem 
and  the  corresponding  solutions  X(y,s,X)  form  a  complete  orthogonal  set 
([7]  Ch.  VII). 

Also,  if  X  4  X  ,  we  have  both 

%(y,s,x)^"(y,s,xn)dy  =  y  %(y,s,x)^(y,s,xn)  pg(y)  -X^]dy 

/%(y,s,X)Z"  (y,s,Xn)dy  -   [%(y,s,X)%^(y,s,Xn)  -X'(y,s,X)^(y,s,xJ] 


X     (y>s,X)%(y,s,Xn)  dy 


and 


00 

o 


-  %'(0,s,X)X(0,s,Xn)  -X(0>s>VX\o,s,\n) 


00 
•■  /  ^(y,s,X)X(y»s,Xn)[s2g(y)  -  xj  dy  . 


Hence 

00 


X'(0,s,X)^(0,s,X  )  -#(0,s,X)%'(0,s,X  ) 

(3.1)    /  ^(y,s,X)%(y,s,Xn)dy  =  

^/o 


X  -  X 

n 


and  so,  if  we  let  X  4  X  ,  we  obtain  by  d'Hopital's  rule 

n 
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SO 


(3.5) 


[%(y,3,xn)]' 


dy  = 


%'H>S>X)    %(o,s,xn)  .  a£(o,s>x)  x»(0js^ 


where, because  of  the  boundary  conditions,  one  of  the  terms  on  the  right-hand 
side  is  zero.  If  we  denote  the  right-hand  side  of  (3»5)  by  co  then  a  complete 
normalized  orthogonal  set  is  formed  by  the  functions 

%(y,s,xj 


(3.6) 


tfn(y,s) 


60 


w 


n 


Now  suppose  that 


(3.7) 

Then 
(3.3) 


CO 


0   (x,y,s)   =  £]  VX,S)  ^n(y,s)   * 

n=l 


So 


I  (x,s) 


n 


^(x,y,s)  Xn(y,s)  dy  . 


Multiply  (3.2)  by  a-  and  integrate  with  respect  to  y  from  0  to  oo»  Then,  since 

n 

K    -  (■Vw-gx, 


and  since  both  a.  and  0     satisfy  the  boundary  conditions  at  y  =  0,  it  follows 
(if  we  assume  that  0  — ^  0  as  y  -^  co)  that 


(3.9) 


ft 


dx 


f  -  (s2f(x)  +0  Yn=  -\(V  5(x-Xo}   * 


Thus,  at  x  -  x  , 
'  o' 


(3.10) 


r   -ix  +  0 
o 


n 


x  -  0 

-*  o 


0,     XT 


r~         —i 

X  +   0 

dl 

0 

n 

dx 

x  -  0 

— ■ ■        — • 

0 

-^JyJ  • 


n'V 
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When  s  is  large,  the  homogeneous  equation  corresponding  to  (3.9) 
has  two  solutions  which  are  asymptotic  to 


K(s2  f(x)  +  X)"l/2  exp  J  £  f   p  f(x)  +  \ll/2 


dx  \    , 


where  K  is  some  constant.  Let  us  denote  these  by  Y  ,  Y   respectively.  We 

n   n 


then  have 

I*  (x,s)  T"  (x  .») 


n 


(3.11) 

Y+  (x  ,s)  Y"(x,s) 


X.  w  (x  <  =0  > 


W  "n  *"o'  '  —  o' 

n 

where  W  =  W  (s)  is  the  Wronskian 
n    n 

dY"        dY+ 

(3.12)  W  (s)   -  Y  -£    -    Y"  -=2   . 
v  n        n   dx     n   dx 

finally,  by  (3.8)  and  (3.6), 

-  00       Y+  (x,s)  r  (xrt,s)%(y,s,\J%(y^,s,Xj 

(3.13)  0  (x,y,s)      =  2Z      -2- 


n       o       'vv''no''n 


n=l  W  (s)  co  (s) 

n  n 


for  x>x.  If  x  <  x  .  then  x  and  x  must  be  interchanged. 
o         o'  o  to 


I4,  Asymptotic  approximations  for  Re(s)  large 

Just  as  in  the  case  of  diffraction  by  a  circular  cylinder,  the 
essential  step  now  is  the  derivation  of  approximation  formulas  for  (3*13) 
that  are  asymptotically  valid  as  Re(s)  ->  co.  For  this  purpose,  the  W.K.B. 
method  could  be  used,  except  that  the  dependence  of  any  X  (n  fixed)  on  s 
is  such  that  the  error  in  the  W.K.B.  approximation  does  not  tend  to  zero  as  s  •*  00. 
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A  more  elaborate  approximation  to  %(y,s,\)  is  therefore  required,  and  can  be 
obtained  by  following  Langer's  procedure^.  For  simplicity,  we  shall  only 
discuss  the  case  of  real  positive  s. 

Let 
(h»l)  X  s  s  (i  , 

If  can  be  shown  that  the  eigenvalues  X. ,  X„,  . ..,  are  in  the  range  u-  >  g(0). 
Hence  we  take  u.  >  g(0).  Then  g(y)  -  u.  changes  sign  when  y  ■  y  «  y  (n),  where 
g(y  )  =  V-t   and  we  can  define  a  real  function  ^  of  y  by  the  equations 


y 

(U.2) 


*  A   Lg^y  *  ■  zi 


2      3/2  A  r     i        -]l/2        ,  * 

3  ^      ■  A   Lg(y  *  •  y      dy  *  y  >  y  > 


■  f     [^  -  g(y 


3/2 /  r    <  q1/2  ,  • 


2   /  ^3/2  /^  ,    t,  H-/'    ,    i  * 


3  (~V       =  /      [i^  -  e(y  )J      d7  >  y  <  y 


• 


M 


Then  ^   has  the  same  sign  as  y  -  y  ,  and 

It  is  easily  shown  that  d^/dy  >  0,  and  that   ^   and  its  derivatives 
up  to  any  required  order  are  finite  at  y  =  y  ,  if  g(y)  is  sufficiently  regular. 
The  behavior  of  ^   as  y  ^oo  depends  on  that  of  yj  for  instance,  if  g  ~ym  for 
large  y  then 

?  ^   y^+2)/3      ^      d?/dy^y(-D/3    # 

Let 


^lAW/3,,     .(g)- 


where  the  Airy  functions  Ai(y)  satisfies  the  differential  equation 
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(I*.*)  Ai     (y)      =     YAi(Y) 

and  is  given  by  (see    [6]) 

.3 


(U.6) 


JOO  o 

Ai(Y)      ■     -    /        cos(  ^-  +  yz) 


dz 


1/2 


Vj 


AN 


As  a  second  solution  of  (li»£)   we  take 


(h.7) 
so  that 

(U.e) 

Now 


ind  so 


Bi(Y)     =      iXl-      ! 


lY3^ 


A1(y)   Hl'(y)  -  ai(y)   bl(y)    -    7    • 


I     .     s^/3  (^  AiCs^  y)      + 


dy\dy/ 


Ai(s2/3  y)    , 


d2X 
dy 


w(^ii-(.^T,tW^)  ^(^) 


-1/2 


Using  (U»3),   (iiolj)   and  (ii«5)   we  can  write  this  as 


ah 
d7 


(U.9) 


[s2    [g(y)   -  |3   +  q(y)"l    X 


q(y) "   UM     ^7  lw] 


As  the  second  solution  of  (U.9)  we  define 


(luio) 


r^ri/2  Bi(s2/^> . 


L^FJ 


-ll^- 
It  then  follows  from  (lu8)  that 

s2/3 
(lull)  X  T  -  X  Y   -   -5-   • 

5,  An  integral  equation  for  % 

We  can  now  derive  an  integral  equation  satisfied  by  X(y)«  In 
the  first  place,  since  both  X  -»•  0,  X  ->  0  as  y  ■»  co,  we  have  from  U.S>) 

*>*  (3.3) 

"X."  X  dy  =  s2  P   (g-n)XX  dy  , 


y  7 

00  _         ^ 


y 

Kence 


/°  X"l  ciy  =   [%'X  -  X  X*]  "  +  /"  [s2(g-u)  +  q)]  X*  dy 


Similarly, 


X*(y)  X(y)  -X(y)  x'(y)  =   /°  x(7l)  q(y1)X(y1)  d^  . 

^y 

y 
X'(y)Y(y)  -X(y)Ay)  -  X'(o)y(o)  -X(o)y'(o)  -  /"Vy^qCy^XCy^dji 


Since  X  can  be  multiplied  by  an  arbitrary  constant,  we  may  suppose  that 

2/3 
XiO)   Y(0)  -X.(0)  y'(0)  =    a—     . 

If  then  follows  at  once  that,  upon  eliminating  X  (y)  £r°m  "the  preceding 
equations,  we  have 

(5.D  X(y)  =X(y)  +^Jy°yX(y)Y(y1)q(y1)X(y1)'ly1  +y^°0X(y1)Y(y)q(y1)X(y1)dy1 

If  we  put 


>  « 
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(5.2)         X*(y)   .  ^[f]1/2,    ,*w   .   [^]-^f^[^] 


-1/2] 


tdy2 


then  (5.1)   assumes  the  form 


2l*(y)      =     Ai(s2/\)    +-§y    (/7Ai(s2/37)Bi(s2/371)q*(y1)X*(y1)dy1 


(5.3) 


s  '        I  ^o 


+     /^AlCe^j)  Bi(s2/3?)   q*(y1)X*(y1)  d7]_ 


>     > 


where  ^,  is  related  to  y.  as  ^  is  related  to  y« 

We  now  make  the  assumption  that  both  the  integrals 

|q*W  U  ,      /    U*|2  ^ 


_oo 


w 


dy 


converge  uniformly  in  p.  in  any  bounded  ii-interval,  and  we  put 

co 

(5.3a)  K     =     K(|i)      -    J         |q*(y)|   dy     . 

-Jo 

o 
Then,  since    |^(y)  |     is  integrable  (according  to  a  fundamental  theorem  of  Weyl) , 

we  have,  by  Schwarz'  inequality, 


00 


kVldy  <  {J***?**         \J^ 


1/2 


l^l2  f     * 


11/2 


<     co 


Hence  an  upper  bound  for  the  integral  terms  in  (5«3)   is 
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n 

7& 


Ai(s2/3^  ) 


max 
0  <  y±<  y 


Bi(s2/3  tl£ 


r7  Iq**,*!^ 

-S  o 


(5.W 


Bi(s2/3^) 


max     Ai(s2/3  ?_)      /"Vx^y.    I     . 

\  >  y  ^y 


2/3  2/3 

Now  for  ^7   >  0,  Ai(s       y  )   is  decreasing,  Bi(s       ^)    is  increasing,  and  their 
product  tends  to  zero.     For   %  <  0,  both  functions  are  bounded.     Thus  there 
exists  a  constant  M,  independent  of  s,(i  and  the  functions  g(u-),  such  that 


Ai(s2/3^  ) 


max 
0  <  y     <  y 


BiCs2/3^) 


<  M  , 


Bi(s2/3<?  ) 


max 

71>y 


Ai(s2/3  XJ 


<    M 


for  all  y  >  Q,    Hence  (5»h)   has  the  upper  bound 

<x> 


(5.5) 


nM 

"273" 


^        |q*X*|dy 


We  can  now  apply  this  estimate  in  two  different  ways.  First,  (5«3) 
and  (5.5)  give 

o 

and  if  this  is  multiplied  by   |q  (y) |   and  integrated  over  (0,oo)   it  follows  that 


*^/* 


q'-ridy    <        (1-^j)"1 


|Ai(s2/3f)    q*|<^     • 
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2/3 

provided  that  s  '  >  nMKj  consequently 

(5.6)  /"  l-ftfl  T    2   «'*    (1-77?)"1    ' 

f  I  2/3  I 

where  M    is  the  upper  bound  of  the  absolute  constant   |Ai(s  '  71  )  |, 
Second,  we  have  from  (5*3)   and  ($•£) 

|X*W    .    Ai(52Ap|    <     jfc  jTU'tfW  , 

and  then  it  follows  finally  from  (5.6)  that,  for  all  y  >  0, 

(5.7)  k*(y>-Ai(s2/3?)|  <  -%«1 


nM  M  K 
8  /J  -  nMK 


6.  Asymptotic  formulas  for  the  eigenvalues 

Let  us  now  take  the  boundary  condition  at  y  =  0  to  be  A.  =  0.  Then 
also  %     =0,  and  so  (5.7)  becomes,  for  y  =  0, 


TtM  m'k 


s  n  -   nMK 


(6.1)  |Ai(s2/3\) 

where,  by  (U»2)  ,  # 

(6.2)  |(-ro)3/2  =  J        t-g(y)]l/2  dy,   |i  -  g(y*)  . 

The  inequality  (6.1)  involves  only  s  and  nj  on  the  right-hand  side 
M  and  M  are  fixed  numbers  and  K  is  a  function  of  p.  which,  since  it  was  assumed 
that  (5o3)  converges  uniformly  in  p.  in  any  bounded  |A-interval,  is  a  continuous 
function  of  11,     It  was  noted  at  the  beginning  of  Section  h   that,  necessarily, 
p.  >  g(0).  Let  us  therefore  suppose  that  g(0)  <  u-  <  g(0)  +  A,  where  A  is  a 
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positive  constant.  Then  if  in  (6.1)  K  is  replaced  by  its  maximum  in  the  num- 
erator and  by  its  minimum  in  the  denominator,  this  inequality  becomes 

nM  m'k. 


Ai(s2/3^o) 


-   775T 


KMK, 


,3  A 


where  K  ,  K?  are  fixed.  If  s  >  (nMK.)  '  )  then  there  exists  a  constant  B 
such  that 

,-2/3 


(6.3) 


Ai(s2/3^o) 


<  B  s 


(8  >  S0) 


As  s  •>  oo,  the  left-hand  side  of  (6.3)  is  0(s         >}  hence  we  must 


2/3 

take  s   V      to  be  close  to  a  zero  of  the  Airy  function  when  s  is  large.  These 

zeros  may  be  denoted  by  -a.,  -a?,...,  where  a.  <  a?  <  •••)  thus 
(6.U)  -?  -  a„  s"2/3  +  0(s"2/3)   . 


a  s 
n 


The  constant  implied  by  the  symbol  0  depends  on  the  order  of  n  of  the  zero, 
but  this  can  be  ignored  if  we  fix  our  attention  on  a  particular  zero. 

Our  final  step  is  new  to  deduce  an  inequality  for  u.  from  (6.10  and 


(6.2) .  These  two  equations  give 
* 


3/2 


&  -  g(y)]1/2  *  -  f  V  *  0(s"2/3»  • 

Now  g(y)  is  by  hypothesis  an  increasing  function  of  yj  hence  we  may  put 
g(y)  -  g(0)  =  Z,   y  =  h(£)  and  write  the  above  equation  as 


(6.5) 


-g(0) 


[n-  g(0)  -S]l/2  h'(S)o% 


o  a 
3  s 


3/2 


+  0(s"2/3) 
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If  we  now  suppose  that 

(6.6)  h'(0)  -  -£—      i      0  , 

g  (o) 

then  (6.£)  is  equivalent  to 


,3/2 
3   s 


|  !s_  ♦  0(3-5/3, 


whence 


H  =  un  -  g(0)  *  [g'(0)]2/3  an  s"2/3  ♦  0(s-l/3) 
or,  by  (li.l) 
(6.7)  \   »  g(0)  s2  +  [g'(0)]2/3  an  sU/3  +  0(s2/3)  . 

This  is  the  required  approximation  for  X  ,  valid  as  s  -^  oo  for  fixed  n. 

Two  remarks  must  be  made  about  this  result.  The  first  is  that  by- 
starting  with  an  iterated  form  of  the  integral  equation  (5.3)  a  better  approxi- 
mation for  X  can  be  obtained.  It  appears  that.  X  can  be  expanded  as  an  asyrap- 

2/3 
totic  series  in  decreasing  powers  of  s   •  A  formal  method  for  obtaining  its 

coefficients  is  presented  in  Appendix  I. 

Secondly,  if  g(y)  -  g(0)  ~  g^  (p  ^  1),  as  y  -f  0,  then  a  formula 

for  X  different  from  (6.7)  is  obtained.  This  point  is  also  discussed  in 
n 

Appendix  I. 

If  the  boundary  condition  at  y  -  0  is  X  =  0,  a  similar  argument 

shows  that  the  n-th  eigenvalue  is  given  by  the  same  formula,  (6.7),  except 

t   i  t 

that  the  a  must  be  replaced  by  the  zeros  a-,   a  ,  ...  of  Ai  (-  a). 
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7<>  Asymptotic  formulas  for  the  eigenf unctions 

We  can  now  derive  asymptotic  formulas  for  the  functions  of  s  which 
appear  in  a  typical  term  of  (3.13).  For  this  purpose  we  write 


(7.1) 


..  n2/3 


s 


■  n  -  g(y*)  -  g(0)  +  [g'(0)]    a   s"2/3  +  0(s2/3) 


withour  specifying  whether  a  is  a  or  a  .  Then  by  reversing  the  argument  of 
the  preceding  section  we  obtain 

(7.2)  7o  -  -  a  s"2/3  +  Ots"**/3)  . 

From  now  on,  we  shall  omit  the  error  terms,  and  write  for  example 

(7.3)  X(y)~(<g)-1/2  Ai(s2/V  . 


We  begin  by  deriving  expressions  for  ^(0),  X  (0)  and  co  .  In  the 
first  place 

so  that 


/-X^ 


a  s 


^7T 


(7.U) 

and  hence 

(7.50 

Also. 


$1    ~    fc'(«] 


L   7Jy=o 


1/3 


X(0)  -  [g'(0)]"V3    Ai(-a)     . 

*  «  ~  &  (%y*  ^i\ , .  1(f)-3/2  i%  Ai(3^  , 


-  a  - 


But,  by  (h.3),   (7.1)   and  (7.2), 


g'(y)  &(y} "  >3 

d7-      =    '2&-g(y)]1/2(-F)1/if  2(-?)3^ 


d2* 


1/2    r 


n  -  g(y) 


1/2 


5(0) 

2^(0]^  a  7& 


L  -f    J 

(i  +  o(3-2/3)]  +  I  W^  [x  +  o(s-2/3)] 


Hence 
(7.6) 


0(1)     . 

*'(0)^32/3    [g'(0)]l/3     Ai'(-a)      . 


Now,  since  X  is  defined  by  (7.1),  we  have  for  example, 

±   X(o,s,x)  =  [g'(o)]2/3  sU/3  ±   X(0,s,X)  +  0(5^/3)  . 


Hence  w  ,  which  is  defined  by  (3.£),  is  given  by 


n 


»n~  -  [gVo^-^Ail^Cg'fo)]-2/3  .+/»  [?(0)]1/3  s2/5AiVn) 
-2/3  [8'(o)J1/3  ii'C^)  fto)]^.^  t'W]^  m'(^) 

and  so,  by  (h.5) , 


(7.7)     «n  ~ 


In  order  to  obtain  % (y,s,X  ) ,  we  require  an  approximation  for  ^ 
when  y  >  y*.  From  (7ol),  y*  is  of  order  s"  '  ,  and  so  for  fixed  y  (y  ^  0) 
and  sufficiently  large  s,  %  is  given  by  (7.3),  where 
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"  /  ^  - 


21*?  -   /  O.N  ^  ^ 


If  we  write,  in  a  notation  slightly  different  from  that  of  Section  6, 
g(y)  =  K,   y  =  h(£),  this  becomes 


3/2 


g(y) 


iVll  =   P  ^_.aV2 


(£  -  H)X/   h  (*)  d£ 

/   ,       27 3/2 
inhere  n  is  given  by  (7ol)  .  Here  — * ,  considered  as  a  function  of  (J.,  is 

continuous  for  u.  =  g(C),  and  its  derivatives  can  be  obtained  by  differentiat- 
ing with  respect  to  [i  under  the  integral  sign,  provided  that  h  (£)  is  suf- 
ficiently regular  and  that  the  resulting  expressions  for  derivatives  of  order 
higher  than  the  first  are  interpreted  as  'finite  parts'  in  the  sense  of 
Hadamard.  Hence, by  Taylor's  theorem, 


2% 


3/2         -j(y)  .  --,.»*.        .     _e(y) 


P-      P      [^Coa^h'WdS-lfg'colJ^as-^Y       i-g(0)]-^h'(,)« 
Jo  Jo 

*   «.-VS)  , 


whence 

3/2 


(7.8)     ilj- 


■/fc,  W  .--  ft',,]'"  »«/  j^  •«.«, 


2/3 

Now  when  s  is  large,  s  '  ^  is  large,  and  so  (7*3)  can  be  replaced  by 

X  (y)  -C^-)"1  2    gnVg  3V6     1A     e^E  3  s  ^V2] 


2"Vi  [g(r>  -  ,]1A  *1/6 


exp 


[-!^3/2] 
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whence,  by  (7*8) 

X(y)  ^      i  />>,-.  ,     1  ,  ,->!  yi,    i  /a      ex?  i-  s^/      [e(y')  -  g(o)J       dy 


5^75  [g(y)  .  g(0)]lA  SV6     «* 


(7.9) 


* ^2/3 ■ sl/3/  w^f) 


Finally,  we  can  obtain  asymptotic  formulas  for  Y"~  (x,s)  and  for 


n 


W  (s)  by  applying  the  W„K«B«  method  to  (3»?)«  We  may  take 


fs2f (x)   +  \J"  exp  J   +     A pf(x')   +  xj         dx 


+ 

r 

n 


whence 


W  (s)  ~   2   . 

n 


Also,  by  (6.7) 


x 


V2  »  ' 

dx    -^  s 


(fU* 


>  ♦  g(o7]1/2  &' 


Hence 


(7.11) 


y;(x,s)y;(xo,s) 

Wn(s) 


'■vy 


4S/  ^wj^4v'^  »^f  rfuS  /j-^} 

2s[f(x)   +  g(0)]VU      [f(xQ)    ♦  gCO)]1^ 
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8.  Approximations  for  0  and  0 


It  was  shown  in  Section  3  that  formally,  for  x  >  x  , 


Y>,s)   l>„,s)  %  (y,s,   O  9C  Cyn,s,Xj 

(8.D     0  =  n  k>  h-  JL— 


00         _ 

n=l 


n    o' 


n  'o'   '  n 


W  (s)  a  (s) 
n  n 


and  that  a  similar  expression  for  0  for  x  <  x  is  obtained  by  interchanging 
x  and  x  •  Suppose  first  that  both  y  >  0  and  y  >  0.  Then  we  can  substitute 
(7.11),   (7.?)   and  (7o7)   into  the  equation  for  0  .     If  we  put 


r 


(8.2) 


(8.3) 


/^Cx')   +  g(0)]l/2  dx'     *Jir\si7)   -  g(0)]l/2  dy' 

So 
ffiiA    ♦  g(0)]-l/2  «.'     -jf  [g(y')    .  6(0)]-l/2  dy' 


(8Ji)  F(l)     -     [fU)   *  g(0)]1/U  ,      G(y)     -     [g(y)  -  g(0)]lA 

£n-     (5'(-aJ*    an[M<-  "n»]2  - 

[g'(o,]2/3 


(8.5) 

then  we  have 


n  8n  enF(x)F(xo)G(y)G(yo)    s2/3 


(8.6) 


exp  |-sr  -  |   [gT(0)]2/3  an6  s1^  ]    [l  +  CXs"1/3)] 
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-1/3 
where  the  error  term  is  0(s  '  )  because  it  is  of  the  same  order  as  the  error 

terms  in  the  expansions  of  the  arguments  of  the  exponentials  in  (7o9)  and  (7.11) • 

It  can  be  shown  that  (8.6)  holds,  not  only  for  real  positive  s,  but  also  for 

all  s  satisfying  Re(s)  >  0. 

The  approximation  (8.6)  cannot  be  substituted  in  (8.1)  indiscriminate- 
ly. In  the  first  place,  the  series  (8.1)  will  not  be  a  solution  of  the  Laplace- 
transformed  problem  unless  it  converges.  It  can  be  shown  to  be  convergent 
(except  when  x  =  x  )  for  real  positive  s.  If  it  also  converges  on  some  line 
parallel  to  the  imaginary  axis  in  the  region  Re(s)  >  0,  and  is  dominated  by 
a  function  which  has  an  inverse  Laplace  transform  defined  by  the  usual  com- 
plex inversion  formula,  then  term-by-term  inversion  is  legitimate.  Now  the 
example  f»0,  g=l+y  which  is  discussed  in  Appendix  II  shows  that  the 
series  may  not  converge  for  Re(s)  >  0  in  the  whole  xy-plane.  But  this  example 
also  suggests  that  term-by-term  evaluation  may  be  permissible  in  a  physically 
important  region,  which  may  be  called  the  'deep  shadow1;  this  region  in- 
cludes a  part  of  the  boundary  y  =  0  which  lies  in  the  shadow.  (The  exact 
definition  of  the  shadow  will  be  given  in  the  next  section.) 

We  now  proceed  on  the  assumption  that  terra-by-term  inversion  is 
permissible.  In  order  to  use  (8.6)  to  derive  an  approximation  for  0  and  then 
one  for  0,   we  must  consider  to  what  extent  the  asymptotic  behavior  of  0  can 
be  inferred  from  that  of  the  0   .  Obviously,  the  factor  e    in  (8.6)  is  a 
•delay  factor',  that  is  to  say  0     is  obtained  by  replacing  the  variable  t  by 
t  -  X    in  the  inverse  Laplace  transform  of  e    0   •  Since  T  >  0,  then 


e    0      must  be  a  Laplace  transform,  and,  a  fortiori;  it  must  tend  to  zero 
as  Re(s)  -^oo©  Hence  (8.6)  can  be  used  only  when 
(8,7)  6  >  0  . 
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It  will  be  shown  in  the  next  section  that  this  inequality  defines  the  shadow, 
so  that  our  approximation  is  suitable  for  the  estimation  of  the  diffraction 
effect. 

Finally,  we  observe  that  the  ratio  of  the  leading  term  of  the  n-th 
mode,  with  n  >  1,  to  the  error  term  of  the  first  mode  is 


^    <  v  «,)  a  > 


0U-1/3) 


which,  because  of  (8.7),  tends  to  zero  as  |s[  -^  00  in  Re(s)  >  0.  In  other 
words,  the  error  term  of  the  first  mode  is  ultimately  greater  than  the  lead- 
ing term  of  any  succeeding  mode*  Therefore  we  can  retain  the  leading  term 
of  the  first  mode  only  when  we  use  (8«6).  We  have  then,  for  y  >  0 

y  >  0,  as  Re(s)  ■*  00 
o 

(8.5)   g  ,      [g'(0)]2/3 l  +  0(fl/3)   exp  f-  ST  -  I^C)]273^  s 

8n  e1F(x)F(xo)G(y)G(yo)     s2/3         ]_ 


If  either  y  -f   0  or  y  ->•  0,  then  the  factors  l/G(y)  or  1/G(y  )  re- 
spectively become  infinite,  and  (8.8)  must  be  modified*  The  requisite  modifi- 
cation is  different  for  the  two  boundary  conditions  at  y  ■  0,  0  =  0  or 
60/dy  =  0.  Suppose  first  that  d0/dy  =  0  at  y  =  0.  Then  0  /  0  if  either  or 
both  y  and  y  are  zero.  Suppose  next  that  y  =  0.  Then  the  term  ^(0,s,X  ) 
in  0^  must  be  replaced  not  by  (7»9)  but  by  (7»5)j  and  (8»8)  must  be  replaced  by 


(8.9)   0  .    l/2,   k'(°}J13 1  +  °[%l/3)e xpJ-sr.|[g'(0j]^6sl/3  [' 

Uti1/2^  Ai(-a1)F(x)F(xo)G(yo)     s1/2 


(y  -  0,  yo  >  0) 
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where  c^  is  the  first  zero  of  Ai  ( -y) .  [The  quantities  f  and  6  retain  the 
meanings  defined  in  (8.2),  (8.3) •]  It  may  be  noted  that  (8.9)  can  be  ob- 
tained from  (8.8)  by  replacing  the  factor  l/G(y)  by  2n1/2  Ai(-  al)  [g'(O)]"1/3  s1/6. 

If  y  >  0,  y  =  0,  then  0  is  given  by  (8.9)  with  y  replaced  by  y. 
If  y  =  0  in  addition  to  y  =  0,  then  0  can  be  obtained  from  (8,9)  by  replac- 
ing l/G(yo)  by  2n"1/2  Ai(_  a^)  fg'(0)]"l/3  s"l/6  : 


(8.10)    0 


2^  F(x)F(x  ) 


1  +  oCs'1^) 

37T^ 


exp  )-   sT 


-  |[g'(o,]2/3  ^6  8V3 


(y  -  y0  -  o)  . 

When  the  boundary  condition  at  y  =  0  is  that  0  =  0,  then  30/dy  can 
be  evaluated  approximately  from  the  leading  term  of  30.. /3y  by  means  of  (7.11), 
(7.9),  (7.7)  and  (7.6):  we  have 


(8.11)    |    - 


_, sM 

l1n1/2Ai  (-o1)F(x)F(xo)G(y( 


.    [>♦  o(i)]^{-^-  |&*(o)]2/\^1/3 


(y  =  o,  yQ  >  0)  , 


where  a.   is  the  first  zero  of  Ai(-y) .  When  y    =  0  with  this  same  boundary 
condition,  we  have  the  trivial  result  0=0,  since  each  0    then  contains  a 
zero  factor.     But  30/dy    retains  a  meaning;  it  is  the  field  due  to  an  'in- 
stantaneous doublet1  at  x  =  x  ,  y  =  0.     The  form  of  an  approximation  for 
this  field  can  be  deduced  by  replacing  y    by  y  in  (8.11).     For  y  =  0  one 


can  then  evaluate  d  0/3y3y  j  it  is  given  by 


(8.12) 


d20 

dydy 


M£    J]  ♦  o(s^)]  ^  L  ST  .  i[g '(o>] 


2F(x)F(x   ) 
o 


1 2/3       .     1/3 
cu.6  s  ' 


(y  =  7Q  =    0) 
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9.  The  geometrical  optics  of  the  problem 

The  inverse  Laplace  transforms  of  the  approximations  (8.8)  -  (8.12) 
can  now  be  written  down  explicitly.  Before  doing  this,  however,  we  propose 
to  consider  the  geometrical  optics  of  the  problem,  as  this  will  enable  us  to 
give  simple  interpretations  of  the  quantities  t",  5  and  the  factors 
[F(x)]-1,  ...,  [g^)]"1. 

It  can  be  shown,  by  a  simple  modification  of  the  proof  given  pre- 

[2] 
viously1-  J  that  the  arrival  time  of  the  disturbance  at  any  point  can  be  de- 
duced from  Fermat's  principle  suitably  modified  when  there  is  diffraction. 
For  the  moment  let  t  =  ^(x,y)   be  the  equation  of  the  front  of  any  distur- 
bance governed  by  our  wave  equation  (1.1).  Then  'f  satisfies  the  eikonal 
equation 

e-»         d?  ♦  (§D2  -  «*>  *  «w  • 

The  rays  associated  with  the  fronts  t  =  f  are  the  characteristics 
of  (9.1),  and  are  defined  by  the  ordinary  differential  equations 

(9.2)  dx   =   £?,    &   .  ?£ 

w'  '  dCT     dx  *    dcr     ay  > 

where  (Ms  a  parameter.  From  (9.2)  we  can  deduce  that 

^2         /»o\2 


dT  dx     dO"  dy     35-  \dxj  \5y/ 

or,  by  (9.1) 

(9.3)  g|     =       f(x)   +  g(y)      . 


? 
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Again,  we  have  for  example 


d2T    dx  a2T     dy 

ax     dcr       axdy   do- 


a2r  ar  +   a2/r   ar 


ax       ax  dx.dy     dy 


i     a    J  fai_ 
2    Sx    i  lax, 


i 


a? 


whence,  by  (9.1) 
(9.U) 


dcr  vax/ 


|f'W 


d_ 

do-  Vdy/ 


=    j  g  (y) 


The  rays  are  also  the  geodesies  associated  with  the  line  element 


(9.5) 


J  fix)   +   g(y) 


2     2 
S.  <{  dx  +  dy 


,2  . 


df   , 


where  the  identification  with  df  is  an  immediate  consequence  of  (9.1),  (9.2) 
and  (9.3),  so  that  f  is  the  geodesic  distance  measured  along  a  ray. 

The  integration  of  (9.2)  and  (9.U)  can  be  effected  by  quadratures. 
First,  we  have 


dT  d  /a' 

dx  cT<rlax"J 


i  »  ,  dx   ar  dt  /art 

2"  f  M   dO"  >       5y  33*  V5y7 


7  s  (y)  a&  » 


so  that 


,ax7 


-  f(x) 


const.,  U-/  -  g(y) 


const. 


Here  ar'/dx,  af/ay  can  be  replaced  by  dx/dO~,  dy/d  0"  respectively,  and  since 
this  can  also  be  done  in  (9.1)  it  follows  at  once  that  we  can  write 


'dx\2 

^dC 


**>♦«*  m 


(9#6)         lasr  ■  fU)  +  ^   (do-r  =  g(y)  ■  ** » 

where  y.  is  an  arbitrary  constant*  Also,  we  can  write  (9.3)  as 


t   -  ^/  [f(x)  +  n]  dO"  +   /  [g(y)  -  tx]  dC> 
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or 

(9.7)  r  =  */[f(x)  +  ^]1/2  dx  +  /[g(y)  -  ^]1/2 


dy 


This  holds  along  any  arc  of  a  ray  on  which  dx/dO",  dy/dG*  are  of  constant  sign, 

and  the  sign  in  (9*7)  must  be  chosen  such  that  both  terms  are  positive. 

Now  consider  the  rays  from  (x  ,  y  ),  the  position  of  the  instantan- 

o   o 

eous  source.  The  equations  of  these  rays  are  obtained  by  integrating  (9.6) 

for  say,  x>x,y<y: 
o      o 

dx _   P°    dy_ 


M)        °"  =  '    [,u',  ♦  $?i  yy    [g(y<, :  ,]iA  • 


X 

o 


(For  x  <  x  ,  the  sign  of  the  first  integral  must  be  reversedj  and 


o 


for  y  >  y  ,  the  sign  of  the  second  integral  must  be  reversed.) 


o 


Now  if  dy/dO"  is  negati  ve  at  x=x,y  =  y,  then  y  decreases  at 
first,  and  there  are  three  possibilities.  If  u-  <  g(0),  the  ray  meets  y  =  0 
at  a  finite  angle.  It  is  then  continued  as  a  reflected  ray.  Defining  the 
•reflection  region*  as  consisting  of  the  points  in  the  plane  which  lie  on  a 
reflected  ray,  we  can  say  that  at  a  point  in  this  region  the  disturbance 
curves  along  a  direct  ray  and  remains  identical  with  the  Hadamard  elementary 
solution  of  (l.l)  until  the  time  t  is  equal  to  the  geodesic  distance  "X 
measured  along  the  reflected  ray  and  the  corresponding  incident  ray.  After 
that,  a  reflected  pulse  must  be  added.  This  pulse  cannot  be  calctilated  from 
the  series  of  modes;  it  can  be  obtained  either  by  converting  the  series  for 
0  into  a  certain  integral  and  evaluating  this  as  an  asymptotic  series  by  the 
method  of  steepest  descent,  or  directly  by  expanding  the  reflected  pulse  as 
a  series  in  half -integral  powers  of  t  -  c  ,  where  T  is  the  arrival  time  of 
the  reflected  pulse 0 
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If,  on  the  other  hand,  u.  >   g(0),  then  the  ray  becomes  parallel  to 
y  «=  0  at  a  distance  y  from  y  =  0  given  by  g(y")  =  p..  Subsequently,  the  ray 
turns  away  from  y  ■  0  again. 

The  third  possibility,  the  transition  case  p.  =  g(0),  gives  the 
glancing  ray.  If  the  shadow  is  defined  as  that  part  of  the  plane  which 
cannot  be  reached  along  a  direct  or  reflected  ray  from  (x  ,y  )  then  the 
shadow  boundary  may  consist  of  the  glancing  ray  beyond  its  point  of  con- 
tact with  y  =  0,  or  it  may  consist  of  a  portion  of  the  glancing  ray  and 
of  the  envelope  of  the  incident  rays0  But  in  any  case  no  reflected  ray 
can  cross  the  glancing  ray* 

Now,  let  (x-,0)  be  the  point  at  which  the  glancing  ray  touches 
y  «=  0.  Then  by  (9«8)  CT  =  6~    where,  since  u.  =  g(0), 


(9.?) 


O; 


dx 


dy_ 


[f(x)  +  g(0)]1/2      /    [g(y)  -  g(0)] 


m 


Let  x  >  ^ ,  and  consider  a  point  (x,y)  in  the  shadow  region  on  the  ray 
touching  y  ■  0  at  x  =  x^,  with  x  >  x_  (see  Figure  1).  Let  0~   be  the 


Figure  1 
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difference  between  the  values  of  C~  at  (x,y)  and  at  (x?,0)  .  Then  clearly 

x  t  J  i 


2 

x. 


[f(X )  +  g(o)]1/2     yo     [g(y )  -  gcoj1/' 


Now  by  (8»3)  we  have,  for  x  >  x  , 

6  =  /*       <*'  /      dy'  /°       dy' 

^    (*.')  *  g(O)]1^  "  /   [,(/>  -  g<0)]^   ^    [g(y')  -  g(0)]^- 
o 

Hence,  by  (9.9)  and  (9.1C), 

(9.11).  6  -   /*   r   ,  te =^7^   . 

^    [f(x  )  ♦  g(O)]1^ 

This  equation  furnishes  the  interpretation  of  6,  and  shows  that  6  >  0  implies 
that  (x,y)  cannot  be  reached  by  a  reflected  ray. 

Next,  we  can  calculate  the  arrival  time  of  the  diffracted  wave  at 
(x,y)  j  it  is  equal  to  the  geodesic  distance  between  (x,y)  and  (x  ,y  )  measured 
along  the  rays  from  (x  ,y  )  to  (x  o)  and  from  (x_,0)  to  (x,y)  and  along  the 
boundary  x_  <  x  <  x  ,  y  ■  0.  Hence  by  (9.7)  and  (9.5),  the  arrival  time  is 

P      [f(x')  +  g(0)]l/2  dx'  +  r°\g{Y)   -  g(0)]l/2  dy  +  /\^)  *   g(0)]l/2dx' 

^x  l/o  \yx_ 

If(x')  +  e(ci)\  dx'  +   /   MvS  -  *(riU  dy'   ■   T  . 


+  P    t(X,)  +  g(0L|l/2  dX'  +   /    [f(y,)  -  S(0)]" 

This  is  the  quantity  T  defined  by  (8.?).  If  x  <  x  ,  similar  considerations 

o 

apply,  and  lead  to  analogous  conclusions. 

Finally,  we  may  observe  that  the  factor  fr(x)F(x  )G(y)G(yOJ~  can 
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also  be  related  to  the  geometry  of  the  rays.  We  take  x  >  x  and  use  the  nota- 

o 

tions  in  (8.2)  and  (8.U),  and  let 

(9.12)  D(x,y)   -  [F(x)G(y)]"1  . 

Then  it  is  easily  verified  that 

This  is  the  transport  equation  associated  with  the  front  t  =  t'Csee  [U],  p  2\\)  . 

From  it  one  can  deduce,  by  an  argument  similar  to  the  one  used  for  the  case 

of  a  homogeneous  medium,  (see  [l]) ,  that  the  product  D  times  the  refractive 

index  |f(x)  +  g(y)J  '  ,  times  the  normal  cross-section  is  constant  along  a 

'tube'  formed  by  adjacent  diffracted  rays.  The  fact  that  D  -^>  co  as  y  -^  0 

is  an  immediate  consequence  of  the  fact  that  y  =  0  is  a  caustic  of  the  diffracted 

rays. 

10.  Approximations  for  0 

The  inverse  Laplace  transforms  of  (8.8)  -  (8.12)  can  be  evaluated 
in  terms  of  the  functions  U,  defined  in  the  author's  paper  on  diffraction  by 

DO 


a  cylinder  UJ.     These  are  defined  by 

c+ico 
(10.1)  Uk(Y,T)     -    ^  /  s"k  exp(s  T  -  Y  sl/3)     ds       . 

^c-ico 


U,  is  zero  for  T  <  0,  and  is  given  for  small  positive  T  by  the  approximation 

rm  *i    n  rv  *t    3(6k'1)A  Y(3-6k)A  T(6k-5) A  pvn  L     ^/2  _\ 

(10.2)     UklY,T)  ~      •   ^1/j     Y        T         exp  I-  -y^   ^yg  I   . 
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Also, 

(10.3)  U2/3(Y,T)  -  32/3  T-V3  H(T)Ai  ^T)"1^  -  ^^   K^  (j^A     H<T> 


and 
(I0.lt) 


3Uk(Y,T) 
5^ 


"U(3k-l)/3  (Y,T)  * 


The  formula  (8.8)  can  be  written,  in  terms  of  (9.12), as 
0   ' 


iL-!L! D(x,y)D(xo,yo)S-2/3  exp  f-sXn-  |[g'(0)]2/3  c^  6s1/3] 


Hence,  from  (10.3) , 

(10.5)    0  ~    *-J£L     D(x,y)D(xo,yo)    HU  K^ 


8n  ^ 


where 
(10.6) 


t  -  r 


i" 


g  (0X0,5) 


3/2 


6(3T) 


TJT—     \   H(T)   » 


is  the  time  counted  from  the  arrival  of  the  diffracted  front*  This  is  valid 
for  small  T,  and  can  be  replaced  for  sufficiently  small  T,  by  the  crude 
approximation 


(10.7)  0 


^[e%)]1/2 


><*,y)D(x0,yo)  ,^-j 


M 


exp 


r 


g  (0)0^6) 


3/2 


6(3T) 


Z&r-    ( 


In  this  formula,  &    and  a.  are  absolute  constants;  6,f  (which  enters  through  T) 

and  D(x,y)  are  functions  of  position  which  can  be  derived  from  geometrical 

t 

optics;  thus  only  g  (0)  appears  as  a  parameter  specifically  connected  with 

the  particular  form,  f (x)  +  g(y)  ,  that  we  have  assumed  for  the  square  of  the 
refractive  index  of  the  medium. 
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The  corresponding  approximations  derived  from  (8.9)  -  (8.12)  are 
as  follows.  If  the  boundary  condition  at  y  =  0  is  that  30/dy  =0,  then 
from  (8.9),  (8.10), 

31/2[g'(0)]l/3     D(xo,yo)        _l/2 

(io-8)  *  -  ,V'J'- 


g,(0)(a'6)3/2 


8n  /£  t^AiC-o^)         F(x) 


^f -—7~i72--  \>    y  =  °>yo>0> 


6(3T) 


(10.9)     0 


^[g'fo)]1^ 


3 

i"/i;  a[  F(x)F (xo)      \ZT/ 


I 


g  (05(^6) 

572 


3/2 


If  on  the  other  hand  0  =  0  at  y  ■  0,  then  (8.11),   (8.12)   give 


(10.10)   || 


[g'(o)] 


U/3 


D(xo^o} 


8n  Ai  (-a^         F(x) 


l£\1/2i 


2~J       TexP<" 


g,(Q)(a16)3/2 


>  » 


y  -  y0  =  0 


7  =  0,  yQ>0, 


(10.11) 


d20 

3ycay 


[g'(o)] 


-1/6 


hn1/2  F(x)F(xQ) 


AVA(3T)^  e*P  <- 


g'(0)(a16)3/2 


6(3T)" 
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Appendix  I.  Asymptotic  Expansion  of  the  Eigenvalues 

The  problem  considered  in  Section  6  was  that  of  determining  the 
behavior,  as  s  -^oo,  of  the  eigenvalues  X  of  the  differential  equation 

(1.1)  A  +   R    S2g(y)]X  =  0 

dy      L        J 

subject  to  the  conditions,  %— >■  0  as  y  ->oo  and  %  *  0  or  d%/dy  ■  0 
at  y  =  0.  The  object  of  this  Appendix  is  to  give  an  alternative  method  of 
attacking  this  problem;  it  is  purely  formal  but  as  such  very  effective. 
Consider  first  the  simple  example 

(1.21         g(y)  =1  +  3^  (p  real  and  positive). 

[it  was  assumed  at  the  outset  that  g(0)  >  6],     Then 

4  *  [x-s2-syu  =  o . 

dy 
Now  let 
(1.3)         I     =  IS«,   a  -  OL-A^I-O  , 

then  we  have,  instead, 

d2*   +  (a  -  f)%     -  0  , 


and  this  has  a  solution  satisfying  the  boundary  conditions  at  y  ■  0  and  at 
V  ■  oo  only  if  a  =  a  when  a_,  o_,  ...  are  a  certain  set  of  numbers  (which 
depends  on  p) .  Hence  we  can  state  at  once  that 

(1.10  X  -  s2+  a  sU/(P+2)   , 

n         n         ' 

that  is  to  say,  we  can  determine  X  explicitly  as  a  function  of  s.  If  p  =  1, 
then  ( J.It)  reduces  to  (6.?)  without  the  error  term. 
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We  can  now  try  to  deal  with  the  general  case  by  a  similar  change 
of  the  independent  variable.   For  simplicity  we  suppose  that 


(1.5)       g(0)  >  0,  g  (0)  >  0,  g(y)  =  gQ+  gly  +  ^y  +  ...  , 

where  e  ■  g   (0)/mJ  and  substitute  (1.5)  formally  in  (I.l).  VJe  also 
assume  an  expansion  for  X: 


(1.6) 


2  °°     -2m/3 
s  )   c  s   ' 

^— »   m 
m=0 


Then  (I.l)  becomes,  in  terms  of  V 
(1.7) 


2/3 


(1.8) 

Now  put 
(1.9) 


d2%      . 


Co  cVl  -  wf +1>  =_V3  U  -  o 


00 


m-0 


-2m/3 


and  substitute  this  in  (1.8).  Then  one  obtains,  by  equating  the  coefficients 

2/3 

of  equal  powers  of  s    to  zero,  a  system  of  recurrence  equations: 


(I.1C) 


(1.11) 


d2Y 


■? 


+  (Cl-g1?)Yo=  0 


d2Y. 


7  +  (V^lPYl=  +(^2?2-  c2>  *c 


s 


d  y  m  • .-, 


-  "W  V) 
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We  can  now  impose  on  each  Y  the  boundary  conditions 

m 

(1.12)  I  ->>  0  as  ^   ->•  °°5   Ym  -  0  (or  dl/d^  -  0)  at  £  -  0  . 

2 
(The  'turning  point'  at  which  g(y)  =  X/s  presents  no  difficulty, 

as  the  substitutions  (1.6)  and  (1.7)  assure  that  its  position  is  essentially 

independent  of  s  for  large  s) . 

A  solution  of  (1*10)  which  tends  to  zero  as  7{  -»  oo  is 

Hence  the  conditions  at  ^  =  0  give 

2/3 

(1.13)  c1   =  a  g^  , 

where  Ai(-a)  ■  0  or  Ai  (-a)  -  0  respectively,  according  to  which  of  the  boundary 

condition  at  71   ■  0  obtains.  Hence 

1/3 
(I.1U)         IQ  =  i/1^    7   -  a)  , 

where  L  is  a  constant.  The  equation  for  Y.,  is  now 
o  1 

(I'15)        ^  *   <•  %/3-  %?>*i  '  V%22  "  °2>  "<^V  -  *    > 


whence 


Tx  -     LjAitg*'3^    -  a)    +  H^Kg*'3^   -  a) 


(1.16)  +     -j^       /    [li(gl/3   ^  »   .  a)Bi(gV3?    -  a)    -  Ai(g^  -  a)Bi(gj/\  *-  a)] 

X     [gs  f  -  c2]  Ai(gJ/3  ^  .  a)    df'      , 
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where  L.,  M.  are  constants.  If  for  example  %  m   0  at  y  ■  0,  then  for  1C  "  0, 

Y  -  M-Bi(-a)   =  0  , 

whence  M.  =  0;  a  similar  argument  holds  if  d£/dy  ■  0  at  y  ■  0,  Thus  M  =  0. 
One  can  also  take  L  *  0,  since  in  the  subsequent  calculation  'L.   /  0  simply 
gives  rise  to  a  series  similar  to  that  which  has  the  factor  L  ,  except  that 
this  one  is  multiplied  by  L^/s  '  •  Then  we  can  take  L  ■  1. 

Now  as  V   -^oo,  the  only  non-vanishing  term  on  the  right-hand  side 
of  (1.16)  is  ^ 

-£j  Bi(Sp,  -  a)  A2  V2-  •$«#*  t  '-  a)]'  d?'. 

This  can  be  written  as 

TIL 


2* 


Bi(gy3- 


Hence  a  necessary  and  sufficient  condition  for  Y.  to  approach  infinity  as 
7  ■»  oo  is  that      m 

f   «%  **"  °2>  (f^Y"  a)]2  «f  ■  0  , 


and  this  determines  c?: 


(1.17)        c 


/%  f  [mU1/3  f  ^  V 


This  argument  can  be  repeated  at  each  stage.  In  fact,  we  can  take 
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^r 


Ym=  -jn 


(1.18) 


h  f\?^    ?'  a)Bi(g^/3  ?  -  a)  -  Ai(g^3  ?  -  ^BKg^^'-  a)] 
1    ^o 


,-^')d?] 


Wm+1-  •*»)  ^^  i- °>  +S(vi?    -  viK 


m-l 


and  we  have  then 


(1.19)     c 


ra+1 


jg  .  _  _  »i  ■*■  . ...  . . 

[iKgj/3  ?  -  aj]2  d? 


where  m  ■  2,3, •• •   • 


It  is  not  clear,  at  this  stage,  when  this  method  is  valid;  nor  have  we  determined 
whether  (1.6)  is  convergent  for  sufficiently  large  s,  or  whether  it  is  an  asymp- 
totic expansion  valid  as  s  -^oo. 
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Appendix  II.  Convergence  of  the  series  for  0  in  a  special  case 

The  simplest  example  of  a  problem  of  the  type  considered  in  this 
paper  is 
(H.1)  f(x)  =0,   g(y)  -  1  +  y  . 

These  functions  satisfy  the  requirements  (1.5)  -  (1.7) •  We  can  then  take,  for 
the  functions  %    defined  in  Section  3, 

X(y,s,\)  =  Ai  fs2/3  y  -  s_1|/3  (X  -  s2)|  . 

Then 

(11.2)  \n  -  s2  +  an  sU/3  , 

where  again  either  Ai(-a  )  =  0  or  Ai  (-a  )  =  0;  this  is  (6.7)  without  the 

error  term.     Thus 

%(y,s,\n)      =     Ai(s2/3y-an) 

and,  by     (3.5), 

0-2/3 

n  n 

where   Cy,  is  defined  by  (8.5).     Also,  we  can  take,  again  in  the  notation  of 
Section  3, 

T*U,s)      -     exp[-x>/r»],     Y;(x,s)      =     expji^/r],     Wn  -     2^     , 

and  so  (3.13)   becomes 

oo       Ai(-a  +  s2/3y)Ai(-a  +  s2/3y  )  f  «  3 A,1/2! 

(11.3)  ?    -C      n-rm ^372 ^    -P  J-     x  -  xo  (s2+  a/A) 

n-1  2e  (s  'J  +  a  )   '  ' 

n  n  l. 

In  order  to  discuss  this  series,  we  must  use  the  asymptotic  behavior 
of  the  Airy  function.  This  is  given  by  the  formulas 
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(II.U) 


Ai^>   -  ^rkjK  exp(-P/2)|i  +  o(Y-3/2)| 


(-n  <  arg  y  <  n) , 


(II.S) 


Ai(r)      =         l/\     l/h     6XP(f    ^)  V  +  0(r"3/27  (  n  <  arg  Y  <  3n), 


where  Y       >  Y       are  to  be  interpreted  as  principal  values.     For  real  positive 


a,  Ai(-a)   is  given  by 


(II.6)  Ai(-a)  =       yf  1/h        <{   sin(|a3/2  +  j)    +  0(<T"-) 


,-3A 


n  '   a 


From  this  one  can  deduce  at  once  that  the  zeros  a  of  Ai(-a)  are  given  by 

n 


(II.7)     f  an/2   =  (n  -  J)*  +  < 


while  those  of  Ai  (-a)  are  given  by 


£>>  %={^-^o(n-^, 


(n  *  1,2,.„.)» 


1 2/3 


(II.8)    |  (03A  ■  (n  -  §»   ♦  0(|) ,   „;  .  (  $  In  -  #  l,"'*  ^ 


(n  *  1,2,...)  . 


Both  these  sets  of  formulas  are  valid  as  n^  coj  the  error  terms  are  numerically 
small  for  n  >  2. 

We  can  now  examine  the  convergence  of  the  series  (II .3) •  For  real 
positive  s,  the  series  converges  for  all  (x,y)  except  on  the  line  x  =  x  , 
since  the  Airy  functions  are  then  bounded.  Next,  suppose  that  0  <  |arg  s|  <  ?  n* 
Then,  for  sufficiently  large  n  (  and  given  x,y,s)  we  have 


*  s2A 


arg  (-a     +     s  '   y) 
n 


<    n 


arg 


<-„♦  ^ 


<    *, 


so  that,  by  (II.U) ,  we  have 
Ai(s2/3  y  -  a  ) 


2^/A     -£3 


£  |  (s2/3 


,3/2l 
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a  similar  expression  holds  for  Ai(s  '  y  -  a  ).  Now  s  '   y  -a  is.  for  laree 

on  n      ' 

n,  in  the  second  or  third  quadrant,  according  as  arg  s  >  0,  arg  s  <  0  respec- 
tively.    Hence 


2   ,  2/3     3/2         2     3/2    l31n/2     L        3     s2/3  y 


2"    +in 
ane 


where  the  upper  sign  applies  when  arg  s  >  0,  the  lower  one  when  arg  s  <  0.  Thus 


Ai(s2/3  y  -  a  ) 


1  (-.   1/2       2/3^ 

1/2  i/h    «p  +ian  y s 


2n  '     a 


n 


n 


Hence  the  general  term  of  (II. 3)  behaves,  as  n  -»  co,  like 

gi-i-S-    «*  £  ^  s2/3[|x  -  x0l  ♦  i(y  *  y0)]j 

Clearly,  a  necessary  and  sufficient  condition  for  convergence  is 

"2/3 


Re  Js 


[|x  -xol  +  i(y  +  yo)^ 


C    >    0     , 


where  the  upper  sign  is  taken  for  arg  s  >  0,  the  lower  one  for  arg  s  <  0. 
Hence  the  series  converges  if 


(II.9) 


y  +  yo  <  |x  -  xj  tan^ 


arg  s 


>   • 


Now  if  (II»3)  is  to  represent  the  Laplace  transform  0  in  the  whole  right-hand 
half -plane  Re(s)  >  0,  then  (II.9)  must  hold  for  |arg  s|  <  j  n,  whence 


x  -  x 


(11.10) 


y  +  y0< 


V? 


Since  also  y  >  0,  the  domains  in  which  the  series  for  0  is  valid  for  Re(s)  >  0 


are  bounded  above  by  two  lines  thr ought  the  image  point  (x  ,  -y  )  inclined  at 


-lit- 


an  angle  of  30°  to  the  positive  and  the  negative  x-axes  respectively,  and, 
below,  by  the  x-axes. 

The  shadow  is  given  by  6  >  0,  where  6  is  defined  by  (8o3)»     In  the 
present  case,  f  =  0,  g  »  1  +  y,  this  becomes 

7*l-  .  /*°aL   >  o  , 


6  .    r  dx'  -  r  sz-  -  r  ^~ 
K*0        J°  v7  J°    V7 


or 


(II.lD  l*-xol     >    2^/7"  +    y^ 


The  relation  between  the  domains  defined  by  (II.9)  and  (11.10)  is 
illustrated  in  Figure  2.  Figure  2a  shows  the  case  y  =  0.  Here  the  'converg- 
ence domain1  (11,10)  contains  the  whole  of  the  boundary  y  »  0  (except  for  the 
point  x  ■  x  )  and,  rather  surprisingly  part  of  the  reflection  region  |x  -  x  |  >  2^p 
as  well  as  a  part  of  the  shadow.  Since  every  diffracted  ray  passes  through 
the  convergence  domain,  it  may  be  possible  in  this  case  to  extend  the  asympto- 
tic formula  for  0  near  the  diffracted  front  by  integration  along  a  ray.  Figure 
2b  shows  a  case  where  y  >  (l  +^/T,)/^/T,     Then  the  lines  (y  +  y  )  -JT   -  |x  -  x  | 
do  not  meet  the  parabolic  shadow  boundary  at  all,  and  so  the  series  expansions 
of  0  can  apparently  not  be  used  in  the  complex  inversion  formula  of  the  Laplace 
transform  in  a  part  of  the  shadow  which  contains  some  of  the  boundary  y  ■  0. 

Clearly,  it  is  necessary  to  obtain  an  alternative  representation  of 
0  which  continues  (II. 3)  analytically  into  the  region 

3    -,   7  ♦  70  1 

i  tan    <  |arg  s|  <  -  n  . 

2         l*-*„i 


-15- 


Formally,  this  is  not  difficult,  as  one  can  replace  (II .3)  by  an  integral 
which  defines  0  for  all  s  with  Re(s)  >  0,  From  this,  an  approximation  to 
$  can  be  obtained  which  is  valid  in  the  region  of  the  xy-plane  in  which 
both  the  incident  and  the  reflected  pulses  occur.  Since  this  turns  out  to 
be  the  geometrical  optics  approximation,  it  can  also  be  obtained  directly. 
In  general,  there  remains  a  region  of  the  xy-plane  (including  at  least  a 
part  of  the  shadow  boundary)  in  which  neither  the  geometrical  optics  ap- 
proximation nor  the  approximation  developed  here  is  valid  (or  useful). 
The  discussion  of  the  behavior  of  0  in  this  region  is  a  problem  which  re- 
mains to  be  solved. 
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